WW Scattering Parameters via Pseudoscalar Phase Shifts 
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Using domain- wall lattice simulations, we study pseudoscalar-pseudoscalar scattering in the maximal 
isospin channel for an SU{3) gauge theory with two and six fermion flavors in the fundamental represen- 
tation. This calculation of the S-wave scattering length is related to the next-to-leading order corrections 
to WW scattering through the low-energy coefficients of the chiral Lagrangian. While two and six flavor 
scattering lengths are similar for a fixed ratio of the pseudoscalar mass to its decay constant, six-flavor 
scattering shows a somewhat less repulsive next- to-leading order interaction than its two-flavor counter- 
part. Estimates are made for the WW scattering parameters and the plausibility of detection is discussed. 



PACS numbers: ll.lO.Hi, 11.15.Ha, 11.25.Hf, 12.60.Nz 

INTRODUCTION 

Technicolor theories based on scaled-up QCD may be 
disfavored due to large contributions to the electroweak S- 
parameter and a chiral condensate too small to account for 
flavor effects fT\. However, walking (T-W\ could account 
for these differences. For a gauge theory with a fixed num- 
ber of colors, walking is believed to set in at a number of 
flavors Nf just below the minimum number required 
to generate an infrared fixed point. Above Nj, the theory 
is conformal in the infrared, and this behavior persists until 
asymptotic freedom is lost. 

Since walking and the onset of infrared conformality in- 
volve strong dynamics, numerical lattice field theory is ide- 
ally suited to study this behavior. To this end, lattice cal- 
culations [5 -9] have been performed with different num- 
bers of flavors, colors, and in different representations, 
leading to evidence for both confining and conformal be- 
havior, with much debate as to the boundaries of these 
phases. The present work builds on two earlier papers of 
the LSD collaboration ifTOl [TTIl , where evidence for both 
chiral condensate enhancement and dynamic reduction of 
the S-parameter were observed for an SU{3) gauge theory 
when Nf is increased from 2 to 6, still well below N^. 

In this work, using domain- wall fermions, we make the 
first lattice connection to WW scattering for these two 
theories. This process, whose longitudinal modes can be 
related to scattering of pseudo-Nambu-Goldstone bosons 
(PNGBs) via the equivalence theorem, is receiving re- 
newed attention. As a direct probe of the physics behind 



electroweak symmetry breaking, WW scattering could be 
an important channel to investigate at the LHC. 

Our approach is to compute the parameters of effective, 
low-energy chiral Lagrangians. We focus on two coeffi- 
cients of the electroweak chiral Lagrangian II21IT3.1 which 
encode the dominant deviations of the longitudinal WW 
scattering amplitude from that of the standard model f\4\ - 
16 1 . These coefficients, in turn, are related to certain coeffi- 
cients of the A^j -flavor hadronic chiral Lagrangian [ TvlfTSll . 
For QCD (Nf = 2), there are two such parameters (li and 
Lattice methods have so far constrained only linear 
combinations of these terms via pion form factors [19 ] or 
the extraction of effective range parameters in 1=2 tttv scat- 
tering iEqI. 

Here we calculate the leading term (scattering length) in 
the effective-range expansion for S-wave, maximal-isospin 
pseudoscalar scattering for the Nf — 2 and Nf — Q the- 
ories. For the linear combination of chiral-Lagrangian pa- 
rameters entering the scattering length, this provides a first 
glimpse of how the reduced running of the coupling associ- 
ated with the increase of Nf affects the dynamics. In sub- 
sequent work, we will compute additional chiral-parameter 
combinations through both D-wave projections and further 
effective range parameters. Together, these will lead to a 
prediction for the two parameters of the electroweak chi- 
ral Lagrangian describing strongly coupled corrections to 
WW scattering. 
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ELECTROWEAK CHIRAL LAGRANGIAN AND WW 
SCATTERING 

The electroweak chiral Lagrangian (EWxL ) allows for 
a systematic description of electroweak scale effects result- 
ing from TeV scale physics ||l2l[n |. This effective chiral 
Lagrangian must respect the SU{2)l x U{1)y gauge sym- 
metry. However, the chiral Lagrangian approximately re- 
spects a larger "custodial" symmetry SU {2)l x SU (2)c, 
which spontaneously breaks to the diagonal subgroup. 
From the EWxL \V2^ lISJ, the dominant contributions to 
WW scattering come from the terms 

+2ig tr((a^Ty. - d,W^)[W^, W,]) 

+a4[tr(V;K)]' + a,[tr{V^V^)]\ (1) 

where = {D^U)U\ U{x) is the unitary matrix field 
that transforms under SU{2)l x SU{2)c (akin to the 
hadronic matrix of Goldstone fields), and the covariant 
derivative of ?7(x) is given by 



(2) 



The and 0^5 terms describe O(g^) corrections to WW 
scattering. 

In Ref. |[2T1 , Eboli et al found that with an integrated 
luminosity of 100 fb~^, and by considering both W 
and Z scattering, the LHC could place 99% CL bounds 
-7.7 X 10-^ < 0^4 < 15 X 10-^ and -12 x 10"^ < 
0^5 < 10 X 10~^. Because the fit to potential LHC data was 
made using only the terms of Eq. ([T]) at tree level, the a- 
parameters are, in effect, defined at a low scale, incorporat- 
ing all radiative corrections including both standard-model 
corrections and new physics. These parameters were also 
constrained by unitarity considerations in Refs. lITSl [T6l . 
Other custodial- symmetry respecting coefficients in the 
EWxL have been constrained experimentally and do not 
lead to appreciable corrections to WW scattering [1151 . 
More recent assessments of LHC constraints on vector bo- 
son scattering were performed in Ref. [22J and Ref. [23|. 

In the limit ^, g' 0, the EWxL reduces to the massless 
two-flavor hadronic chiral Lagrangian [17 J, as illustrated in 
Fig.[T] where 



«5 = j + Oig) 



(3) 



where i\ and ^2 are the next-to-leading-order (NLO) 
Gasser-Leutwyler coefficients of the hadronic chiral La- 
grangian with only derivative couplings. 
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FIG. 1: 

Comparison between limits of the hadronic chiral Lagrangian 
with Nf flavors 1 17 , 18 1 and the electroweak chiral Lagrangian 
Id Hi I . Symbolic quantities are defined in the text. 



To exhibit the flavor-dependent dynamics giving rise to 
0^4 and 0^5, an extension to the multi-flavor hadronic chi- 
ral Lagrangian, with non-zero fermion masses is appropri- 
ate. Theories with multiple flavors lead to additional Gold- 
stone degrees of freedom, uneaten by the W and Z bosons, 
which, with a small mass, become PNGBs. Their contribu- 
tion to physical phenomena are parametrized by the many 
low-energy constants (LECs) of the multi-flavor chiral La- 
grangian. The presence of a finite fermion mass is also 
essential for the lattice simulations employed here. 

For theories with Nf > 4 massive fermions, there are 
9 LECs in the NLO hadronic chiral Lagrangian, denoted 
by Lo-8 mi- The LECs L4_8 multiply terms proportional 
to the fermion mass, while L0-3 multiply terms that are 
independent of the fermion mass. One way to relate these 
LECs to the EWxL is by assigning electroweak quantum 
numbers to one fermion doublet among the Nf fermions, 
leaving the others neutral. 

Although our lattice simulations are carried out with a 
common mass m for all the fermions, it is helpful for our 
discussion to temporarily assign a mass to the fermion 
doublet with electroweak quantum numbers and a mass 
^ to the remaining electroweak singlets. The for- 
mer must be taken to zero, while the latter may or may 
not be taken to zero. The Nf — 1 PNGBs are then sep- 
arated into the 3 (with mass denoted by M^d) which be- 
come massless when rrid and are eaten by the gauge 
bosons, the 4:(Nf — 2) composed of one electroweak- 
doublet fermion and one electroweak- singlet fermion (with 
mass denoted by M^s), and the (Nf — 2)^ composed of 
only electroweak-singlet fermions (with mass denoted by 
M55). All except the 3 eaten modes can also get masses 
from standard-model and beyond- standard-model interac- 
tions not included here. 

The relation of the 9 NLO LECs of the general-flavor 
hadronic Lagrangian to and 0^5 in the electroweak La- 
grangian is mapped out in Fig. [T] "High-energy physics" 
including the PNGB particles composed of one or two of 
the electroweak-singlet fermions of mass are integrated 
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out to describe observed low-energy physics on the order 
of the W and Z masses. Specifically, it is conventional to 
start with a set of 9 LECs L\[ii) defined at some scale ji 
small compared to the breakdown scale of the chiral La- 
grangian. They are independent of PNGB masses, being 
defined in the zero-mass theory, and infrared finite because 
of the presence of the scale ji. 

Now suppose that the two-flavor Lagrangian of Eq. (4) 
refers to the two electroweak-doublet fermions. The two 
coefficients, /^l ^2' ^^^^ t)e defined at some scale 
/i, and in the limit M^^ = 0. Since the PNGBs of mass 
Mfis and Mss have been "integrated out" (integrated into 
and I2) these two coefficients do depend on these masses. 
At NLO and with the PNGB masses small compared to 
the breakdown scale, they are related to the L\{ii) by the 
addition of chiral logarithms arising from a loop of PNGBs 
composed of one EW-doublet fermion and one EW-singlet 
fermion PT8l[24ll: 



ds) 



-2Ll{^i) + ^UM + 2Ll{^i) 



24(327r2 



log 



Ml 



+ 



Nf 



12(327r2 



log 



Ml 



(5) 



The coefficients of the the chiral logarithms are propor- 
tional to the expected counting factor of 4(A^j — 2). 

These quantities are related to the corresponding 
electroweak LECs al{ii,Mds) and a^di^Mds) by 
Eq. ([3]). Then, a set of /^-independent electroweak 
LECs a^{Mdd,Mds) and a4{Mdd, M^s) can be defined 
by adding chiral logarithms arising at NLO from a loop 
of PNGBs composed of electroweak-doublet fermions of 
mass m^: 



a,{Mdd,Mds) 



ds) 



1 



ds) 



384^2 

1 

192^ 



log 



ML 



log #,(6) 

where we have dropped 0(g) corrections. These elec- 
troweak LECs have incorporated into them all strong- 
interaction effects including the PNGBs made purely of 
EW-doublet fermions, given at NLO by the explicit chiral- 
logarithm terms. 

Comparison with the and employed by Eboli et al 
II2TII requires the inclusion of standard-model corrections 
not considered here. Then the above log(Mj^) dependence 
can be separated and incorporated into these corrections 
(allowing the limit Mdd to be taken there) since the 
resultant NGBs are eaten by the W and Z. We therefore 
subtract from Eq. ^ the one-loop contributions arising 
from PNGBs of mass Mdd together with a Higgs boson 
with reference mass Mh- From the first of Eq. ([6]), we 
subtract {-l/3847v^)[log{MlJMfj) + 0(1)], and from 



the second {-l/1927v^)[log{ML/Mfj) + 0(1)], where 
the 0(1) constants are determined by the aforementioned 
one-loop calculation. The remaining quantities are then 

_ lli^.Mds) [log^ + 0(l)] 



a^{MH,Md 



a^{MH,Md 



Ul^.Mds) [log^ + 0(l)] 



3847r2 

Ml 



4 1927r2 

(7) 

The subtracted LECs a^iMn, Mds) and a^^Mn, Mds) 
are directly analogous to the conventionally defined S pa- 
rameter, as employed in Ref. [1 1 J. They can be determined 
experimentally and compared to lattice computations if the 
additional, standard-model corrections, with the same Mh, 
are also taken into account. For Nf > 2, they remain sen- 
sitive to the PNGB masses Mds, and therefore the limit 
rris ^ can be taken only if additional physics provides a 
mass for these modes. 



LATTICE FRAMEWORK 

For the lattice calculation being performed here, all Nf 
fermions are assigned a common mass m, and, as a result, 
have a common PNGB mass Mp. The goal is to compute 
the resultant L\{ii) through pseudoscalar scattering, and 
then make use of Eqs. 5-7 with Mp — Mds — Mdd, to 
determine the quantities a4{MH^ Mp) and a^(MH, Mp) 
as a function of Mp. For Nf = 2, the limit Mp ^ is fi- 
nite. For Nf > 2, the residual log Mp dependence means 
that the limit m ^ can be taken only if interactions not 
included here generate a mass for these modes. 



MAXIMAL-ISOSPIN PSEUDOSCALAR SCATTERING 

The study of pion scattering is well established in ef- 
fective field theory, where the leading order (LO) result 
was first calculated by Weinberg [ 25 1 and the NLO chiral- 
perturbation-theory result was calculated by Gasser and 
Leutwyler [ 18 1. Our focus here is on the scattering of two 
identical pseudoscalars within an SU (2) l x SU(2)p sub- 
group of the flavor symmetry, taken here to refer to the 
fermion doublet assigned electroweak charges. This pro- 
cess is often referred to as "1=2" scattering in the case of 
two flavors. Here, we refer to it as "maximal-isospin" scat- 
tering for Nf fermion flavors with degenerate masses. 

Low-energy pseudoscalar scattering can be parameter- 
ized in terms of the phase shift, 6. On the lattice, the ac- 
cessible quantity is cot 5, where |A:| is the magnitude 
of the pseudoscalar 3-momenta. The S-wave projection of 
this combination has a convenient effective-range expan- 
sion for small momenta given by 



\k \ cot 6 



1 Mlrpp[\k 



+ 



app 



Ml 



+ (9 



\kr 

Ml^ 



(8) 
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where app is the scattering length (in the particle physics 
convention), rpp is the effective range, and Mp is the 
pseudoscalar mass. The expansion is valid for \k\ <C Mp. 

The NLO chiral expansion for the scattering length in 
the MI channel takes the form Il26ll27l 



Mpttpp — 



M2 



167rF2 



1 + 



M2 



167r2F2 



2(iV/ - 1; 



+ ^(iV/)log(— 



where = 2m(^^)/F^, with F and (ipip) the pseu- 
doscalar decay constant and chiral condensate at m = 0, 
and where 



(2 



Upp 



-25Qn'[{Nf-2){Ll{^,)-Ll{^^)) 
+LS(m) + 2LU/x) + 2L^(m) + L^(m)]. 



Each L\{ii) is defined at a scale in the MS scheme. The 
quantity bppifJ^) + A{Nf) log(M^//i^) is /i-independent. 

Similar chiral expansions of the pion mass and pion de- 
cay constant are given by 1271 



Ml, 



M^i 1 + 



M^ 



167r2F2 



1 

^m(m) + irr log 



Nf 



11- 



M^ 



167r2F2 



¥p{^i) 



Nf , M2 



where 



bp{p) 



+2L;iii) - mil)] 
647r'[NjLl{f,) + Ll{ii)]. 



(12) 



The expressions for Mp and Fp are also //-independent. 

Only the LECs Lq-L^ contribute to WW scattering, 
(Eq. ([5])), since only they multiply operators that survive 
as m ^ 0. But in general the scattering length computed 
here will give us the combination of parameters appearing 
in bpp (Eq. ([TO])), which include L4 - Lg. More data is 
required to separate these sets for general Nf. 

For Nf — 2, however, 

5^p(m) = -2567r2 [LKp) + 2L^(/x) 
+2mii) + Ll{fi)] 
^-128Tr^[r,{li) + i;ili)]. (13) 

Both and l^ili) are now Mp-independent since 

the only PNGBs are the three composed of EW-doublet 
fermions, and therefore the chiral-logarithm terms in Eq. 



^ vanish. From Eq. Q, we then have for the Nf 
case 

1 



aiiMn) + a^iMn] 



5127r2 



[log 



Vpp{ll) 

+ 0(1)] 



128.2 • '''^ 
Thus, for Nf — 2, di lattice calculation of Mpapp as a 
function of M, with a fit using Eq. ([9]), can directly deter- 
mine a/^^Mn) + 55 (Mj:/). 
(9) For Nf = 2, it is also worth noting that if the analysis is 
taken to NNLO order in jF^, and Ss can be deter- 
mined separately, employing only the quantities app, Mp 
and Fp. The NLO LECs introduced above multiply single 
chiral-logarithm terms at this order. While promising, the 
primary difficulty is having enough data within the chiral 
regime to constrain these terms accurately. 



LATTICE DETAILS 



(^10) Lattice calculations are performed using 32^ x 64 
domain-wall lattices with the Iwasaki improved gauge ac- 
tion with a fifth dimensional length of = 16 and a 
domain- wall height of mo — 1.8 HOl [TT]|. By using 
domain-wall fermions, the calculation preserves exact fla- 
vor symmetry and chiral-breaking lattice spacing artifacts 
are Lg suppressed (rUres <^ This is performed for 
TV/ = 2 at /3 = 2.70 and TV/ = 6 at /3 = 2.10, where the 
beta values were chosen to match the IR scales of both the- 

Iories and correspond to an inverse lattice spacing roughly 
five times the vector meson mass. For both Nf — 2 and 
Nf — 6, five mass points are analyzed with m — rrif + 
(11) rUres. where ruf = 0.010, 0.015, 0.020, 0.025, 0.030. 



FINITE VOLUME METHOD 

In Euclidean space, scattering phase shifts can be ex- 
tracted on the lattice only by calculating the total energy 
of two hadrons in a finite volume [28|. In practice, this 
is accomplished by calculating the four-point correlation 
function and exploring the long time behavior 



tr(7r(t)^)-tr(7r(t))^ 



^ e 



-Eppt 



(15) 



where the zero-momentum projected operator Ti{t) is re- 
lated to the pseudoscalar two-point function by Cp{t) — 
tr(7r(t)), and Elp = 4(|fcp + Mj,), with \k\ being the 
center-of-mass scattering momentum. In this work, we 
restrict ourselves to S-wave scattering by projecting each 
pseudoscalar correlator onto zero momentum. For a finite 
volume, only discretized values of the scattering interac- 
tion are allowed. This scattering momentum is related to 
the phase shift by [ 28 1 



\k\ cot 5 



4^2 



S\ 



(16) 
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where the function S{r]) is given by the regularized zeta 
function [291 



(17) 



CALCULATION AND FITTING 

The calculation yields two quantities: the four-point 
pseudoscalar correlation function, Cpp(t), and the two- 
point correlation function Cp{t). The Fp values used in 
this work were previously calculated in Ref. [ 10 |. The long 
time behavior of the two correlators differ due to the extra 
backward propagating pseudoscalar that exists in Cppii) 
and as a result, the long time behavior of the midpoint- 
symmetric correlators is given by 



Cpp{t) A + Bcosh{Eppt) 
Cp{t) Ccosh(Mpt). 



(18) 



For the pseudoscalar mass, a hyperbolic cosine fit was per- 
formed, where 



cosh(Mp) 



Cp(t + 1) + Cp(t-1) 



2Cp{t) 

A different fitting form was used for Epp, namely 
Cpp(t + 2) -Cpp{t-2) 



(19) 



cosh(£^ 



ppj 



2{Cpp(t + l)-Cpp(t-l))' 



(20) 



Constant fits were performed to the inverse hyperbolic co- 
sine of Eq. ([19]) and Eq. ( [2Q| ) over a window of roughly 
40 time slices in the plateau region. In addition to statisti- 
cal errors, systematic effects due to placement of the fitting 
window were also examined by varying the window by ±2 
time slices on each side. 



SCATTERING RESULTS 

In our lattice calculations, |A:|cot5 is extracted from 
Eq. ([16]). For Eq. ([9]) to be applicable, the effective range 
expansion, Eq. ([8]), is needed and \k\'^ /Mp must be suf- 
ficiently small. The data suggest that the approximation 
Mp/\k\ cot 6 ^ MpGpp is valid within 9%, and thus 
Eq. ([9]) should be applicable to the lattice results. Our 
results for pseudoscalar scattering in the MI channel are 
given as a function of m in the bottom graph of Fig. [2] 
We first analyze the Nf — 2 results and then discuss the 
difference between Nf — 2 and Nf — Q. 
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FIG. 2: Plots of M|,/2m, Fp and Mp/m\k\coi5 
Mpapp/m versus m (in lattice units). The error bars are sta- 
tistical plus systematic. The red circles represent the Nf — 2 
data and the blue squares represent the Nf — Q data. The fits for 
A^j = 2 are made using only the solid red points. 



Nf = 2 



To determine bpp{fi) for Nf = 2, we carry out a com- 
bined fit using NLO chiral perturbation to Mpapp (Eq. 
([9])), M|, (Eq. ([11])), Fp (Eq. ([TT])), and the chiral con- 
densate. The latter three quantities were computed and an- 
alyzed in Ref. [10], where it was concluded that for the 
range of m values employed there, NLO chiral perturba- 
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tion theory provides an acceptable fit allowing a reliable 
extrapolation to m = 0, determingjhe extrapolated decay 
constant F and chiral condensate (^^). 

TABLE I: Results (in lattice units) of a combined NLO fit for 
Nf = 2 to Mpapp, Mp, Fp, and the chiral condensate. 





m/=0.01-0.015 


mj =0.01-0.02 


m/=0.015-0.02 


Ref. 1 10 




80(10) 


110.7(9.0) 


198(60) 


91(29) 




5.20(28) 


5.38(15) 


5.09(58) 


5.70(27) 




-2.36(21) 


-1.74(24) 


0.5(2.2) 


-2.22(62) 




-17.32(88) 


-16.89(59) 


-13.3(3.8) 




F 


0.0220(16) 


0.0229(10) 


0.0262(35) 


0.0209(41) 




1.049(74) 


0.885(39) 


0.65(10) 


0.99(17) 


xVdof 


16/2 


83/6 


13/2 


26/4 



For the present fit, we choose ji — F c:^ 250 GeV (In 
Ref. fTOl , the scale ji — AtvF was used.) The (combined) 
fits for Ml /2m, Fp and Mp/m\k\ cot 5 2^ Mpapp/m 
are shown in Fig. [2] Only the points ruf — 0.01 — 0.02 
(shown in solid red) are used in the fit. For all quantities, 
NLO chiral perturbation theory again provides an accept- 
able fit as in Ref. [10|. 

The fit parameters are shown in bold face in the central 
column of Table [l] In addition, fit parameters for two other 
m f ranges are shown, and used to estimate systematic er- 
rors. From the table, the quantity of interest here is 

b'ppil^^F) = -U87T'[l[{t,^F) + lUii^F)] 



-16.89 ±0.59tj^j^ 



dof 



83 



(21) 



where we have used the central value from the fit range 
rUf = 0.01 — 0.02. The errors are statistical plus system- 
atic. This result is consistent with previous lattice simula- 
tions of 1=2 pion-pion scattering Il20[|30t^ . other 2-1-1 lat- 
tice QCD determinations of these LECs ElllOl, and QCD 
phenomenology. Lighter mass ensembles will be required 
for future calculations to achieve higher precision. 
From Eq. (14\ , we then have, for Nf = 2, 



(3.34 ±0.17+;^;°^ 



71 J 



10" 



[log# + 0(l)] 



128.^ ' '''^ 
where the errors are again statistical plus systematic. Re- 
call that the expression in the second line arises from hav- 
ing removed the one-loop contribution arising from the 
eaten Goldstone bosons and a Higgs boson with reference 
mass For a wide range of M^, it is O(10~^) or 

smaller. 

Once the additional one-loop standard-model correc- 
tions with the same Mjj are included, the above result can 
be compared directly with LHC data. If the standard-model 
corrections are also of order 10~^ or smaller, then the full 



set of 0(^^) contributions will fall comfortably within the 
Eboli et al bounds based on 100 fb~^ of LHC data. Fu- 
ture measurements will have to be more precise to compare 
meaningfully with Eq. ([22]). 



Nf = 6 

The Nf = Q data for Mp/m\k\ cot 6 2^ Mpapp/m 
displayed in Fig. [2] lies statistically on top of the Nf — 2 
data. A NLO fit using Eq. ([9]) would nevertheless lead to 
a different value for h^pp^ji) because the chiral-logarithm 
term in this expression has a much larger coefficient, grow- 
ing linearly with Nf. However, the larger coefficient indi- 
cates that it is very unlikely that Eq. ([9]) can be employed 
in the Nf — Q case for the existing range of m values. 
This point was already made in Ref. [10] where Fp, Mp 
and the chiral condensate were computed. The NLO ex- 
pression for Fp (Eq. ( pT] )) also has a chiral-logarithm term 
with a coefficient growing linearly with Nf. 

Nevertheless, our results for Nf — Q do indicate an in- 
teresting trend, which can be seen more easily by plotting 
the data in a different way. If the NLO chiral expansion Eq. 
([9]) is reorganized in terms of the physical values Mp and 
Fp, it takes the form 



Mpapp 



Ml 
IGttFI 



N] 



1 + 



Ml 



IGn^Fl 



^A\Nf)\ogy 



fMl 



(23) 



where 



2{l -Nf + Nj) 
Nj 

-2567r2 [Lliii) + 2Ll{fi) + 2Ll{p) 



+Ll{f,)-2Ll{i^) 



Ll{^l) + 2Ll{^l) 

(24) 



With Mpapp expressed in this way, the coefficient of the 
logarithmic term remains finite as Np increases, and the 
LEC combination h'pp{ii) contains no explicit Nf depen- 
dence. 

The above expansion, while not a priori more reliable for 
Nf = 6 than Eq. ([9]), suggests that a plot of Mpapp ^ 

Mp/\k \ cot 6 versus the physical quantity Ml/Fl could 
be revealing. This is done in Fig. |3j where the solid-color 
points correspond to the range mj = 0.01 — 0.02. A small 
upward shift of the Nf = 6 points relative to Nf = 2 is 
indicated. Since a negative value for app corresponds to 
a repulsive interaction, the data indicate that the Nf = Q 
theory is somewhat less repulsive than the Nf — 2 theory 
for pseudoscalar scattering in the maximal-isospin channel. 



7 




(Mp/Fp) 



FIG. 3: Plotof Mp/|A;|cot^ ^ Mpapp vs. {Mp/Fpf. The 
error bars are statistical plus systematic. The red circles represent 
the two-flavor data and the blue squares represent the six-flavor 
data. The dashed line is the LO xPT result (zero parameter fit). 
Larger negative results correspond to more repulsive scattering. 



The dashed line, representing the LO expression 
— Mp/167r^Fp, is a reasonably good first approximation 
to the data for both iV/ = 2 and Nf = 6. For Nf = 2, the 
data show that the effect of the NLO term is to make the 
interaction more repulsive. The quantity in square brackets 
in Eq. ( [23] ) is positive and of order unity within the range 
shown. A fit to just Mpapp with ji — F leads to the 
value h'^p{ii = F) = -4.67 zb 0.65tJ;^^. Clearly there 
is some cancelation between this term and the chiral loga- 
rithm. Nonetheless 



this h'pp value (when combined with 
the 6^ and 6^ values in Table |l]) is consistent with the 6pp 



value in Eq. ([21 

For Nf — 6, the data is even closer to the LO dashed 
line, suggesting that NLO perturbation theory in the form 
of Eq. [23] might again be reliable, 
used to fit the Nf 



If this expression is 
6 data, then the quantity in square 
brackets is again positive and of order unity within the 
range shown, but somewhat smaller in magnitude than for 
Nf — 2. Since we don't yet know the precise value of F 



in lattice units for Nf - 
the scale ji — 0.023a" 
h'^p{lJi = 0.023a-^ - 



= 6, we carry out the NLO fit using 
1 (F for Nf = 2). The fit leads to 
F) = -7.81 ± QAQtlil larger 
in magnitude than for Nf — 2. There is now more cance- 
lation between this term and the chiral logarithm than for 
Nf = 2. 

The above values of h'pp emerge from a fit of Eq. ( [23] ) 
to each of the three lightest data points (corresponding to 
ruf 0.01 - 0.02), with a fixed choice 0.023a-^ 
F. A plot of the resultant value of h'pp versus m (Fig. |4]), 
shows that h'pp{ii = 0.023a"^ ^ F) is relatively inde- 
pendent of m for both Nf — 2 and A^j = 6 as expected 
if NLO perturbation theory is reliable. The evident shift 
going from A^j = 2 to A^j = 6 is interesting since this 
quantity is contains EEC's that enter into WW scattering 
through Eq. ([24]). 



On -4 
(N 
(N 
O 
O 
II 

^ -6 



-10, 



• N =2 

■ N,=6 







i 



!^ W'W W'tt w tt'W wii W'T:^ w iiiw 



'F 



'1 



0.005 0.01 0.015 0.02 0.025 0.03 



FIG. 4: Chiral parameter h'pp versus fermion mass m for Nf 
2 and Nf 6. 



It is not yet clear whether this fit can be trusted for 
Nf — 6, but even if it can, the resultant value for 
hpp{iJ. = 0.023a~^ 2^ F) determines only the combi- 
nation of LECs in Eq. ( [24] ), which includes L\{ji) values 
not directly relevant to WW scattering. Further calcula- 
tions will be necessary to isolate a4(MH^ Mp = M^s) 
and as{MH,Mp = M^s) (Eq. 0). This will then de- 
scribe the effect of beyond-standard-model physics for a 
range of PNGB masses Mp. 



SUMMARY AND DISCUSSION 

Using lattice simulations, we have computed 
pseudoscalar-pseudoscalar scattering in the maximal 
isospin channel for an SU (3) gauge theory with two and 
six fermion flavors in the fundamental representation. 
Our calculation of the S-wave scattering length was then 
related to the next-to-leading order (NLO) corrections 
to WW scattering through the low-energy coefficients 
of the chiral Lagrangian. For Nf = 2, our result for 
the scattering length agreed with previous calculations, 
showing an increase in repulsion due to the NLO correc- 
tions. For WW scattering, we obtained an estimate for 
^4{Mh) + a^{MH) (Eq. ( [22] )) describing deviations from 
the standard model. 

Six-flavor scattering showed a somewhat less repulsive 
NLO interaction than its two-flavor counterpart for a fixed 
ratio of the pseudoscalar mass to its decay constant. The 
range of fermion masses employed so far does not allow a 
clearly reliable use of chiral perturbation theory. Also, the 
appearance of more terms in the hadronic chiral lagrangian 
for six flavors does not allow the extraction of only the 
combination of parameters entering WW scattering. Fur- 
ther simulations of additional low-energy scattering param- 
eters at lower fermion-mass values will be required to com- 
plete this study. 
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